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Abstract 

For an investigation of the physical properties of gravitational fields 
the observation of massive test particles and light is very useful. The 
characteristic features of a given space-time may be decoded by study¬ 
ing the complete set of all possible geodesic motions. Such a thorough 
analysis can be accomplished most effectively by using analytical meth¬ 
ods to solve the geodesic equation. In this contribution, the use of elliptic 
functions and their generalizations for solving the geodesic equation in 
a wide range of well known space-times, which are part of the general 
Plebanski-Demianski family of solutions, will be presented. In addition, 
the definition and calculation of observable effects like the perihelion shift 
will be presented and further applications of the presented methods will 
be outlined. 


1 Introduction 

The observation of massive particles and light is a very important tool for ex¬ 
ploring the features of gravitational fields and also for tests of general relativity. 
The motion of massive and massless test-particles is described by the geodesic 
equation, which is a coupled system of ordinary differential equations depen¬ 
dent on the metric of the considered gravitational field. A wide range of ex¬ 
actly known solutions of Einstein’s field equations possesses certain symmetries, 
which allow to decouple the geodesic equations. Here we discuss metrics within 
the Plebanski-Demianski family of solutions (see [I]), which is a seven param¬ 
eter solution with mass, rotation, cosmological constant, electric and magnetic 
charge, NUT charge, and acceleration and which comprises the Schwarzschild 
and the Kerr metric as special cases. In this family of solutions the equations 
of motion considerably simplify due to the separability of the Hamilton-Jacobi 
equation (for lightlike geodesics and, if the acceleration vanishes, also for time¬ 
like geodesics). 
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Due to this simplification, the analytic solutions for the complete set of 
geodesics in Schwarzschild space-time was already found in 1931 by Hagihara 
[2] in terms of elliptic functions. With essentially the same methods also the 
geodesic equations in Kerr-Newman-Taub-NUT space-times and subcases can 
be handled (and also an additional acceleration for massless particles). For 
the case of the Kerr metric, this was first done in the equatorial plane (see [3] 
for a review) and later, after the introduction of Mino time [4] which allows 
to fully decouple the equations of motion, by Fujita and Hikida for bound or¬ 
bital motion [5], For a nonvanishing cosmological constant, the structure of 
the equations of motions is more complex but can still be solved analytically 
by using hyperelliptic functions as demonstrated for the Schwarzschild-de Sit¬ 
ter metric [6] and for general axially symmetric space-times [7j. Here, we will 
explain these general methods to analytically solve the geodesic equations in 
Kerr-Newman-Taub-NUT-de Sitter space-times (and the C-metric for lightlike 
geodesics). 

For observational purposes also explicit expressions for the deviations of 
relativistic orbits from the Kepler orbits are of interest. Here we concentrate on 
the observables for bound orbital motion, namely the periastron shift and the 
Lense-Thirring effect. However, in a strong gravitational field concepts like the 
orbital plane and the orbital ellipse are are no longer valid. A fully relativistic 
treatment of these effects in the Kerr gravitational field was given by Schmidt [Sj 
and combined with the Mino time by Drasco and Hughes [Sj as well as Fujita and 
Hikida (5j. In this article, we will show how these concepts can be generalized 
to the above mentioned space-times. 


2 Equations of motion 


The motion of test-particles is described by the geodesic equation 


d 2 x p 


dx p dx a 


+ = o 

ds 2 pa ds ds 
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where T%, = \g pa {d pg> 


+ d a g pa — d a g pa ) is the Christoffel symbol and fi = 
0,1, 2, 3. This system of coupled ODE’s can be simplified if the underlying 
space-time has certain symmetries. In particular, for the Plebanski-Demiahski 
solutions (with vanishing acceleration for massive test-particles), there exist 
four constants of motion which can be used for decoupling: the normalization 
constant e = liras' witlu e = 0 for light and e = 1 for massive test-particles, 
the energy E, the angular momentum L in direction of the symmetry axis, and 
the Carter constant K. 

For the considered family of solutions in standard Boyer-Lindquist coordi¬ 
nates the radial and latitudinal equations of motion can be reduced to the form 

m 



P(x-,p) 


( 2 ) 
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where x is the radius or the (squared) cosine of the latitude, A is an affine 
parameter, the Mino time [4], P is a polynomial in x of degree 2g + 1 or 2g + 2, 
0 < i < g is an integer, and p = {p±,... ,p n } is a set of parameters of the 
space-time and the test-particle. For example, in Kerr space-time we get in 
geometrized units (G = 1, c = 1) [31 H] 

^ = ((r 2 + a 2 )E — aL) 2 — (r 2 + a 2 — 2r)(er 2 + K) =: R(r) , (3) 

(10 =^K 1 -0(K-ea 2 Z) + (aE(l-Z)-L) 2 ], (4) 

where £ = cos 2 6 and all quantities are normalized such that they are dimen¬ 
sionless. 



3 Algebro-geometric methods 


The equation of motion © should be solved for x(X), i.e. we consider the in¬ 
version problem 


f x x l dx 
Jx o y/Pixyp) 


(5) 


where a;(Ao) = Xq are initial values. The solution of this problem should be 
independent from the chosen integration path. This implies that for a closed 

path 7 with u := <f f dx ^ 0 the solution x(X) has to have the period w, 
7 V p (x;p) 

= A — Ao — ui. This can be taken into account automatically if the 

x\p) 

equation © is considered as an algebraic curve 

w 2 =P(x-,p), ( 6 ) 


r x x % c 
Jx° 


where w = x' 1 ' 4?. For the considered space-times one of the following two 
situations occur 

• P is of order 3 or 4: then © is an elliptic curve of genus 1, 

• P is of order 5 or 6 : then © is an hyperelliptic curve of genus 2. 

Topologically, (hyper-)elliptic curves can be considered as Riemann surfaces. 
The genus g corresponds to the number of ’holes’ in the Riemann surface, see 
Figure [I] If there are g holes this implies that there are 2 g independent closed 
integration paths whose integrals do not vanish and, therefore, the solution s(A) 
needs to have 2 g independent periods [TO]. 

In the case of an elliptic curve © can be reduced to the Weierstrass form 
by a rational substitution, 

w 2 = 4i 3 — g 2 x — g 3 . (7) 
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Figure 1: Riemann surface of a genus one curve (left) and a genus two curve 
(right) with 2 or 4 independent paths. 


In this standard form, w and x are parametrized by the Weierstrass elliptic 
function and its derivative, w = p'(z) and x = p{z). In the above example of 
Kerr space-time, the substitution for the radial equation of motion for a timelike 
geodesic is r = + r K , where R(r K ) = 0 and aj = fijjyr ^ (4 -p ( r ic)- The 

resulting equation is given by 0 with 




1 (1 2 3 2 

93 — Yg I -aia 2 a3 — — a 2 — a 0 a 3 


( 8 ) 


The analytical solution for the timelike radial equation in Kerr space-time is 
then given by 


r(A) 


a 3 

4p(A - c) 


02 
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+ r K , 


(9) 


where c = c(ro, Ao) is a constant which depends only on the initial conditions. 

With a complete analogous procedure all equations of motions which re¬ 
duce to elliptic curves may be solved. This includes geodesics in Schwarzschild, 
Reissner-Nordstrom, Kerr-Newman, and Taub-NUT space-times as well as the 
C-metric (see e.g. [IT]) for massless particles. 

The situation gets more complicated if the equation of motion is described by 
a hyperelliptic curve of genus two. This is due to the fact that the solution has 
to have four independent periods, which is impossible for functions in a single 
complex variable. Therefore, is is necessary to consider a function in g = 2 
complex variables. However, as we have only one degree of freedom, we need 
to restrict the space of variables again to a one-dimensional submanifold. This 
is done by considering the equation of motion as part of the Jacobi inversion 
problem for g = 2, which is to solve the system of equations 


W\ 


w 2 




dx 
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for x\, X2 as functions of w\, W2- If P is transformed to a standard form 
P(x) = 4x 5 + X]i=o diX 1 the solution to this problem is known in terms of gen¬ 
eralized p-functions, X1X2 = @12(1x1, u’2), x\ + X2 = @22(1x1,1x2) [ID]- Here 
pij(w 1,1x2) = — logcr(wi, W2) with the generalized a function. By let¬ 

ting X2 go to infinity, we can get rid of the second integral on the right hand 
side and simultaneously restrict (ix\,W 2 ) to the one-dimensional sigma divisor, 
i.e. the set of zeros of the two-dimensional cx-function. The solution for x = X\ 
is then given by 

x = — — (wi,W2) , with a(wi,W 2 )= 0 , (11) 

where <r, denotes derivative of <r w.r.t. the i-tlr variable. For example, in the 
case of Schwarzschild-de Sitter space-time, the solution for the radius r in terms 
of the angle ip is given by |6j 

r(<p) = -M—(f(tp),ip), with a(f(tp),p)=0. (12) 

01 

This solution method can be applied to all geodesic equations which reduce 
to hyperelliptic curves. This comprises the radial and latitudinal equations of 
motion in Kerr-de Sitter space-time as well as all the de Sitter-versions of Kerr- 
Newman-Taub-NUT space-times and subcases [7|. For lightlike geodesics also 
the C-metric with an additional cosmological constant can be treated this way. 


4 Observables 


For the bound orbital motion of massive test-particles there are two main gravi¬ 
tational effects: the periastron shift Ap and the Lense-Thirring effect Alt- The 
first is a precession of the orbital ellipse within its orbital plane and the latter a 
precession of the orbital plane itselfbJ In this picture, the periastron shift is de¬ 
fined as the angle between two consecutive periastrons and the Lense-Thirring 
effect as the angle between two consecutive minimal latitudes. Mathematically 
this means that they are given by the difference between the period of r(p) or 
9(p), respectively, and 27r. If 2 u r denotes the period of r(p) and 2 cog the period 
of 9(ip) then 

Ap = 2w r — 2 sign(L) 7T , Alt = 2we — 2 sign(L) tt , (13) 


where the sign of L is included to distinguish between prograde and retrograde 
motion. 

In Schwarzschild space-time, where r(p) is directly known, the exact ana¬ 
lytical expression for the periastron shift is given by 


4 LK(k) 

V ( E2 - - n) 


-2t r, 


(14) 


the original paper na a combined effect was computed: the timely variation of the 
ascending node and the argument of periapsis. With ’Lense-Thirring effect’ we refer only to 
the first correction. 
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(15) 


where K(k) is the complete elliptic integral of the first kind, 

dx 


K(k ) = I" 
Jo 


1 0 y/(l-x 2 )(l-k 2 x*) ’ 
which is implemented in standard mathematical software. Here k 2 = ri y a ~ r p| 

r p\ r & ^ 1 ) 

with the zeros 0 < ry < r p < r a of R(r) = • 

The expression for the perihelion shift can be generalized to axially sym¬ 
metric space-times, where only r(A) and </?(A) is known. Writing <p(A) as a part 
linear in A, given in form of an infinite Mino time average T^, plus oscillatory 
deviations mu, we may use A(y>) = T^ip. Then the period of r(ip) is given by 


r(X(ip + 2A r T ip )) = r(T V + 2 A r ) = r(T V) = r(\(<p)), 


(16) 


where 2A r is the period of r(A). Accordingly, 2uj r = 2A r T v and the perihelion 
shift can be written as 


Ap = 2A r T ip — 2 sign(T) n . (17) 

In a totally analogous way the Lense-Thirring effect can be found. If 2 Ag 
denotes the period of 6( A) we may write again 9(\(<p + 2AgT v )) = 9(\(<p)) and, 
therefore, 


A lt = 2A g T tp - 2 sign(L) tt. 


(18) 


In the case of Kerr-Newman-Taub-NUT-de Sitter space-times, A r g and T v 
are given in terms of (hyper-)elliptic integrals. For the elliptic case, they can 
be rewritten in terms of the three standard Jacobian elliptic integrals A'(fc), 
E(k), and n(n, k). For the case of genus two or higher, to our knowledge such 
a standard form does not exist. A possible generalization of the first Jacobian 
elliptic integral would be 



i A it l ldt 


(19) 


where A is a vector of length g, which reflects the fact that there are g indepen¬ 
dent differentials of the first kind, and k is of length 2g — 1 . 

For example, in terms of these integrals the perihelion shift for Schwarzschild- 
de Sitter space-time is given by 

C4 R i (r4 —r3) A (^1 , k "2 , ^3) 

A P = _ i3, ' ; ' - - 2 -, (20) 

sjL 2 Kr 3 (ri - r 0 )(r 4 - r 2 )(r 5 - r 4 ) 


where ro < 0 < r\ < r 2 < < r 4 < r§ are the zeros of the defining polynomial 

with 7"3 = 7' m in and r 4 = r max for small positive A and 


2 _ro(r 4 -r 3 ) , 2 _r 2 (r 4 -r 3 ) 

1 r 3 (r4-ro)’ 2 r 3 (r 4 - r 2 ) ’ 


H = 


^ 5(^4 ~ r 3 ) 
r3<,r 5 - r 4 ) ' 


( 21 ) 
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5 Outlook 


The methods presented here are powerful tools for the analytical integration of 
the geodesic equation in a wide range of space-times. Beside the space-times we 
focused on here, geodesics in higher-dimensional spherically symmetric space¬ 
times mm and the Meyers-Perry space-time m can be treated. It may also 
be possible to extend these methods to the equations of motion in space-times 
with given multipole moments like the Erez-Rosen space-time. 

The hyperelliptic curves which we used here to represent the equations of 
motions are a special case of abelian curves, which allow for higher orders of w in 
(l6j) and mixed terms. These more general curves appear for example in Horava- 
Lifshitz and Gauss-Bonnet gravity, which may be represented by quartic curves 
of the form (w — P{x)) 2 = Q{x). A generalization of the presented methods to 
these cases is in preparation. 

Analogously to the analytic expressions presented here for observables of 
bound orbital motion also the bending of light and the gravitational time delay 
may be considered. Linked to that, we plan the development of an analytical 
timing formula for pulsars orbiting a black hole. 

Concerning the numerical calculation of the analytical expressions, the com¬ 
plete elliptic integrals can very efficiently be computed by using the arithmetic 
geometric mean. This can be generalized to genus two hyperelliptic integrals, 
see [16’, 17 r 
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